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Abstract 

We describe a method for determining the dispersion of magnetic field vectors about 
large-scale fields in turbulent molecular clouds. The method is designed to avoid inac- 
curate estimates of magnetohydrodynamic or turbulent dispersion - and help avoiding 
inaccurate estimates of field strengths - due to large-scale, non-turbulent field struc- 
ture when using the well-known method of Chandrasekhar and Fermi. Our method 
also provides accurate, independent estimates of the turbulent to large-scale magnetic 
field strength ratio. We discuss applications to the molecular clouds OMC-1, M17, and 
DR21(Main). 

Subject headings: ISM: clouds — ISM: magnetic fields — polarization — turbulence 



Introduction 



Chandrasekhar k, Fermil (119531) use d the dispersion of starlight polarization vectors about con- 



tours of Galactic latitude (IHiltner 
velocity dispersion 



19491 ) 



together with estimates of gas density and line-of-sight 
to determine the strength of the magnetic field in the arms of the Galaxy. The 
same technique, "The Chandrasekhar- Fermi, 'CF', method", has been applied, with modifications, 



ature, wavelengths, sensitivities, and resolutions (e.g., 


Lai et al. 


2001, 


2002, 


2003; 


Crutcher et al. 


2004; 


Houde 


2004; 


Girart et al. 


2006; 


Curran & Chrvsostomou 2007 


)■ 



The basis for deriving field strengths from dispersion measurements is the same for observations 
of Galactic arms or molecular clouds: in either case dispersion decreases as the field strengthens. 
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But in the case of the Galactic arms, the dispersion is due to magnetohydrodynamic (MHD) waves; 
the displacements are perpendicular to the direction of propagation. In the case of turbulent dis- 
persion in molecular clouds, there is no preferred direction. The turbulent component can be in any 
orientation. 

Moreover, in dense clouds, the field may have structure due to effects such as differential 
rotation, gravitational collapse, or expanding H II regions; i.e., structure not accounted for by the 
basic CF analysis. Consequently, dispersion measured about mean fields, assumed straight, may be 
much larger than shoul d be attributed to MHD waves or turbulence. Disp ersion measured about 
model large-scale fields (|Schleuning)ll998l : lLai et al.ll2002l : iGirart et al.ll2006l ) that give approximate 
fits to a polarization map will result in better estimates but still give inaccurate values of the 
turbulent component, since they are unlikely to perfectly match the true morphology of the large- 
scale magnetic field. In this paper we describe a method for determining magnetic field dispersion 
about local structured fields, without assuming any model for the large-scale field. This method 
also provides accurate, independent estimates of the turbulent to large-scale magnetic field strength 
ratio. 

We begin (Sj2]) by discussing difficulties one must overcome in order to infer turbulent structure 
from polarization maps, regardless of large-scale effects. In $31 we present the method and in $H 
we give applications to the molecular clouds OMC-1, M17, and DR21(Main). Detailed derivations 
resulting in the relations and functions used in the aforementioned sections, as well as the data 
analysis, will be found in the appendices at the end of the paper. 



2. Difficulties in Deriving the Turbulent Structure from Polarized Emission 



Turbulent velocities of gas motion within and between clumps of mat erial along the line of sight 



can often b e inferred from the widths and centers of molecular lines (e.g.. iKleiner h Dickmanlll984l . 



19851 . 1 19871 ) . But dust polarization measurements of dispersion in magnetic field direction do not 
separate contributions from either volume elements located along the line of sight or across the area 
subtended by the telescope bea m. Hence the measured angular dispersion tends to be a smoothed 
version of the true dispersion ([Myers fc Goodmanl Il99ll : IWiebe fc Watsonl 120041 ) . Nonetheless, a 
corresponding average of the dispersion remains and is measurable; for a given object observations 
will thus reveal a higher degree of dispersions when they are realized at an accordingly higher spatial 
resolution. 

A potentially fruitful line of attack for estimating magn etic field strengths relies on comparisons 
of observed and simulated map s of the net polarization (e.g.. lOstriker et al i l200ll : lHeitsch et al.l[200ll 



Falceta-Goncalves et al.ll2008l ). If the simulations are computed for the resolution, column density, 
and other characteristics of the cloud under study, and if they are computed for several models of the 
key variables (e.g., field strength and turbulent fraction), then one can find the model gi ving the best 



fit to the observations. A valid simulation must also take into account temperatures (jVaillancourt 



3 



2002 ) and grain alignment efficiencies in different environments (IHoang Lazarianl 120081 ) . The 
comparisons are facilitated if both the observations and the simulations are presented in tables of 
Stokes parameters, so that each can be analyzed in the same way. The various modifications of the 



CF method that have been used to relate net dispersion to field strength ( e.g., lOstriker et al.ll2001 



Padoan et al 



2001 



Heitsch et al 



2001 



Kudoh fc Basu 



2003 



Houdell2004l ) are, in effect, first-order 



substitutes for simulations. 

But a meaningful comparison between simulations and observations can only be achieved if a re- 
liable estimate of the spatially averaged angular dispersion can be secured experimentally. It would 
therefore be advantageous if a more general method, which does not depend on any assumption 
concerning the morphology of the large-scale magnetic field, were devised. The method we describe 
in the following section allows for the evaluation of the plane-of-the-sky turbulent angular dispersion 
in molecular clouds while avoiding inaccurate estimates of the turbulence and corresponding inac- 
curate estimates of field strengths due to distortions in polarization position angles by large-scale 
non-turbulent effects. This method can lead to valid estimates of magnetic field strengths only un- 
der conditions such that the Chandrasekhar-Fermi method can be properly applied: a smooth, low 
noise, polarization map, precise measured densities and gas velocities that are moderately uniform, 
and an adequate accounting of the integration process implicit to polarization measurements. This 
latter aspect will be addressed in a subsequent paper. 



3. A Function to Describe Dispersion about Large-scale Fields 

Consider a map precisely showing the angle $(x) of the (two-dimensional) plane-of-the-sky 
projected magnetic field vector B (x) at many points in a molecular cloud. We obtain a measure of 
the difference in angle, A<£ (£) = $ (x) — $ (x + i), between the N(£) pairs of vectors separated by 
displacements £, also restricted to the plane-of-the-sky, through the following function 



1/2 



_g[$(x)-S(x + *)] 2 j 



1/2 



where (• • •) denotes an average and t = \t\. The square of equation (ED) is also often referred to 



as a "structure function" (of the second order in this case; see iFalceta-Goncalves et al.ll2008l : iFrisch 
19951 ). but for our applications we shall refer to it as the "dispersion function" and assume that it 
is isotropic (i.e., it only depends on the magnitude of the displacement, £, and not its orientation). 
We seek to determine how this quantity varies as a function of I. 

To do so, we will assume that the magnetic field B (x) is composed of a large-scale, structured 
field, Bo (x), and a turbulent (or random) component, B t (x), which are statistically independent. 
We also limit ourselves to cases where S < £ <C d, where 5 is the correlation length characterizing 
B t (x) and d is the typical length scale for variations in Bq (x) . 
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Focusing on Bo (x) we would expect its contribution to the dispersion function to increase 
(since (A$ 2 (£)} is positive definite) almost linearly starting at £ = and for small displacements 
£ -C d, as would be expected from the Taylor expansion of any smoothly varying quantity. We 
denote by m the slope characterizing this linear behavior. We also expect a contribution from 
the turbulent component of the magnetic field B t (x). This contribution will vary from zero as 
£ — > (when the two magnetic field vectors are co-aligned) to a maximum average value when the 
displacement exceeds the correlation length S characterizing B t (x) . More precisely, we expect that 
the turbulent contribution to the angular dispersion will be a constant, which we denote by 6, as 
long as £ > 5. These two contributions must be combined quadratically, since the large-scale and 
turbulent fields are statistically independent, to yield 



(A$ 2 {£)) ~ b 2 + m 2 £ 2 (2) 

when 5 < £ <i 

A more formal and rigorous derivation of equation (j2j| is established in Appendix [A] under the 
further assumptions of homogeneity and isotropy in the magnetic field strength over space. Although 
these assumptions are unlikely to be realized across molecular clouds, this level of idealization is 
necessary to allow us to gain insights on, and some quantitative measure of, the importance of the 
turbulent component of the magnetic field in molecular clouds. 

In reality, the measured dispersion function from a polarization map will also include a con- 
tribution, <7m (£), due to measurement uncertainties on the polarization angles 3> (x) that must be 
added (quadratically) to equation J2|). The square of the total measured dispersion function then 
becomes 



(A^ 2 (£)) tot ^b 2 + m 2 £ 2 + a 2 M (£). (3) 

when 5 < £ < d. The function (A$ 2 (^)} tot , not (A$ 2 (£)), is the one calculated from a polarization 
map (from an averaging process similar to equation [T]), and will thus contain separate components 
due to the large-scale structure (i.e., m£), the turbulent dispersion about the large-scale field (i.e., 
b, the quantity we wish to measure), and measurement uncertainties (i.e., cm (£))■ 

If there were no turbulence and no measurement uncertainties, then, for f < (i the measured 

2X 1 / 2 



dispersion function would be a straight line with zero intercept, ( A<3? {£) \ = m£ (see Figure [TJ 

\ / tot 

Curve A). Taking the measurement uncertainty, % (^), into account, the line would be displaced 
upward as specified by equation (J3]) (Curve B, where cjm was assumed to be independent of £). 
Likewise when we next consider turbulence, the curve will again be displaced upward in the same 
manner (Curve C) except at values of £ below the angular resolution scale at which the observa- 
tions were made (Curve D), or below the turbulent correlation scale 5 (Curv e E). Theoretical and 



observational estimates of 5 for molecular clouds are on the order of 1 mpc (ILazarian et al.l 12004 : 
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Li fe Houddl2008l . respectively), well below the size of the telescope beam with which the observa- 
tions presented in this paper were obtained. Although it has not yet been feasible to resolve 8, it is 
now feasible to determine the turbulent dispersion at scales comparable to the approximately linear 

portion of ( A$ (£) 



2 y/ 2 

\ / tot 



Notice that o"m {£) can be accurately determined through the uncertainties on the measured 
polarization angles of each pair of points used in the calculation of <^A<I>(£) 2 ^ , and by then 

subtracting its square to obtain (A&(£) 2 \. As the number and precision of the vectors improve, 



equation j2]) can be fitted to the data for 8 < £ <C d, and the intercept at I = provides us with 
the turbulent contribution, b 2 , to the square of the angular dispersion. 

The Chandrasekhar-Fermi method for evaluating strength of the plane-of-the-sky component 



of the large-scale magnetic field (jChandrasekhar k, Fermilll953l ) implies that 



* * vT> (4) 

where 8B stands for the variation in the magnetic field about the large-scale field Bq, a (v) is the 
one-dimensional velocity dispersion of the gas (of mass density p) coupled to the magnetic field, and 

is the Alfven speed. It is further assumed that the dispersion, <r$, in the polarization angles $ (x) 
across a map is given by 

SB 

a$ ~ — . (6) 

-DO 

The combination of equations j5|), and j6|) allows for the aforementioned determination of the 
plane-of-the-sky component of the large-scale magnetic field strength as a function of p, a (v) (de- 
termined from the width of appropriate spectral line profiles), and <t<j> (determined from polarization 
measurements). 

It is shown with equation (IA24j) in Appendix [A] that the ratio of the turbulent to large-scale 
magnetic field strength is given by 



1 /2 

It is therefore apparent that we should make the correspondence (B 2 ^ — ► 8B and that 



- 6 - 



B Q ~ V(2-6 2 )4tt /3 ^ 
o-(v) 



~ a/87T/5 



where the last equation applies when B t <^ Bq. The fact that the turbulent dispersion, b, is to be 
divided by approximately v2 before being inserted the Chandrasekhar-Fermi equation is readily 
understood by the fact that (neglecting the contribution of the large-scale field) 



(A&(£)) = ([$(x)-<IKx + £)] 2 ) 
= 2 (<0> 2 > - (*) 2 ) 
= 2cr|, 

when £ > 5. Since we also know that (A<1? 2 (^)) = 6 2 at these scales, we then find that b 2 = 2a 2 > , 
which is consistent with equations j6j) and J7|). 

It should be noted that the combination of equations JTJ and j8]) allows, in principle, for the 
determination of both the large-scale and turbulent magnetic fields' strength from polarization and 
spectroscopy data. 



4. Applications to the Molecular Clouds OMC-1, M17, and DR21(Main) 



Using data from the polarimeter Hertz (Dowell et al. 1998) at the Caltech Submillimeter 
Observatory at 350 fj,m, we have measured dispersion functions for the mole c ular clouds OMC- 
1, M17, and DR21(Main). Thes e data are di scussed in detail in Irlourie et~ai1 (ho04h for OMC-1, 
Houde et all (l2002h for M17, and iKirbvl (bofld ) for DR21(Main). Figure [2] shows the results for all 
sources. More details on the data analysis will be found in Appendix iBl 

1 /2 

For each object, we show (A<3? 2 (^)) over the cloud along with the best fit from equation 
J2]) using the first three data points to ensure that t <C d, as much as possible. The measurement 
uncertainties were removed prior to operating the fits to the corresponding data sets. The turbulent 
contribution to the total angular dispersion is determined by the zero intercept of the fit to the 
data at £ = 0. The net turbulent component, 6, is 0.18 ± 0.01 rad (10.4 ± 0.6 deg), 0.12 ± 0.02 rad 
(6.8 ±1.3 deg), and 0.15 ±0.01 rad (8.3 ±0.3 deg) for M17, DR21(Main), and OMC-1, respectively. 

Although large variations in density within the observed regions prevent a reliable estimate in 
the field strength at precise locations, it is still possible to give some average value for the large-scale 
and turbulent field strengths. To do so we use representative line width measurements from H 13 CO + 
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J = 3 — > 2 detections withi n the three clouds. For OMC-1 and M17 we have used the corresponding 
measurements published in lHoude et al.l (|2000i ) (more precisely, an average of the variances obtained 
at the two positions listed for M17), while for DR21(Main) we have used previously unpublished 
data. This molecular species is well suited for this as the effective dens ity needed for line detection 
with the aforementioned transition (n e g ~ 10 5 cm -3 , see lEvanslll999n is close to the densities at 
which dust continuum emission is detected at the measured wavelength. Also, the corresponding 
spectral lines are likely to be optically thin (like the dust continuum) and an ion molecule such as 
this one is better coupled to the magnetic field (and the dust) than corresponding neutral species 
(e.g., H 13 CN for th e same rotational transition) over the whole turbulent energy density spectrum 
(|Li fc Houddl2008l ), Therefore, using a density of 10 5 cm -3 and a mean molecular weight of 2.3 we 
obtain the results shown in Table U As a simple comparison, the values of dispersion shown in 
the table are approximately three times lower than would be obtained if one naively calculated the 
dispersions about the global mean field (i.e., the field direction defined by the mean of all polarization 
vectors in the corresponding map). More precisely, we get dispersions of 27.2, 21.0, and 26.8 degrees 
about the global mean field orientation for M17, DR21(Main), and OMC-1, respectively. 

We wish to emphasize the fact that the quoted values for Bq could not be precise to better 
than a factor of a few due to a lack of precise gas density numbers. Moreover, the values for the 
large-scale magnetic field strength we derived are up to an order of magni t ude h igher than those 
obtained with other observational means (cf., the results of ICrutcher et al.l (|1999l ) for OMC-1 and 
M17 using CN Zeeman measurements). These high values are in part the result of the smaller 
angular dispersions obtained using our technique as compared to more common methods used 
when applying the Chandrasekhar- Fermi equation (e.g., model fits to large-scale fields). One must 
keep in mind, however, that the process of signal integration through the thickness of the cloud 
and across the telescope beam that is inherent to polarization measurements has not been taken 
into account. We will show in a subsequent publication h ow this situation is rectif i ed when these 
considerations (and others) a re carefully taken into account ( Myers fc Goodmanlll991 : Ostriker et al 



2001 : Wiebe Sz Watsonll2004 ). Nevertheless, the turbulent to large-scale magnetic field strength ratio 
is precisely evaluated through our equation (JT]). 



5. Summary 

We have described a method to estimate plane-of-the-sky turbulent dispersion in molecular 
clouds while avoiding inaccurate estimates of the turbulence and corresponding inaccurate estimates 
of field strengths due to distortions in polarization position angles by large-scale non-turbulent ef- 
fects. The method does not depend on any model of the large-scale field. We plot a "dispersion 
function", the mean absolute difference in angle between pairs of vectors as a function of their dis- 
placement £ and show that this function increases approximately linearly for displacements greater 
than the instrument resolution, greater than the correlation length, 5, and less than the typi- 
cal length scale, d, for variations in the large-scale magnetic field (§4). We emphasize that this 
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method can lead to valid estimates of magnetic field strengths only under conditions such that the 
Chandrasekhar-Fermi method can be properly applied: a smooth, low noise, polarization map, pre- 
cise measured densities and gas velocities that are moderately uniform, and an adequate accounting 
of the integration process implicit to polarization measurements. This method, however, provides 
accurate estimates of the turbulent to large-scale magnetic field strength ratio. 

Although the resolution of the instruments now available are not adequate to directly determine 
the correlation length, S, one can still determine the dispersion in the fields at scales where 5 < £ <C d 
for the angular dispersion function. We have successfully done this for the OMC-1, M17, and 
DR21(Main) molecular clouds. 

We thank Shantanu Basu for helpful discussions. This work has been supported in part by 
NSF grants AST 05-05230, AST 02-41356, and AST 05-05124. L.K. acknowledges support from 
the Department of Astronomy and Astrophysics of the University of Chicago. M.H.'s research 
is funded through the NSERC Discovery Grant, Canada Research Chair, Canada Foundation for 
Innovation, Ontario Innovation Trust, and Western's Academic Development Fund programs. J.E.V. 
acknowledges support from the CSO, which is funded through NSF AST 05-40882. 



Let us define the total magnetic field B (x) as being composed of a deterministic, Bo (x), and 
a turbulent (or random), B t (x), components such that 



A. Dispersion Relation Derivation 



A.l. Analysis in Three Dimensions 



B(x) =B (x)+B t (x). 



(Al) 



These quantities have the following averages at points x and y 



(Bo (x)> 
(B (x)-Bo (y)> 

(Bt (x)> 
(B (x) • B t (y)> 



B (x) 

B (x) • B (y) 







(B (x))-(B t (y)) = 0. 



(A2) 



We will further assume homogeneity in the field strength over space. That is 



<B§(x)> 
<B?(x)> 



{K(y)) = Bl 

(B t 2 (y)> = (B t 2 >. 



(A3) 
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Let us now consider the quantity 



(COS [A$ 3 D (£)]) = 



(B (x) -B(x + l)) 

[(B* (x)>(B2 ( x + /))] 



1/2- 



(A4) 



The quantity A3>3D {£) is the angle difference between two magnetic field (or polarization) vectors 
separated by a distance £, the average of its square is the function that we wish to evaluate through 
polarization measurements (albeit in two dimensions, see §A,2j) . Using equations (lAlj) and (|A2j) we 
find that th e numerator of equation (|A4j) (i.e., the autocorrelation of the total magnetic field; see 
Frischlll995l ) becomes 



(B(x)-B(x + ^)> = £ 2 + (B (x) 



Ml 

£-(e,.VfB (x) 



.71=1 



+ (B t (x)-B t (x + £)>, (A5) 



where we used the Taylor expansion 



EMI 
-(e,-V) n B (x), 



71=1 



(A6) 



with ei the unit vector in the direction of £. 

If we introduce d the scale length characterizing (large-scale) variations in Bo and we consider 
situations where £ = \£\ <C d, then we would expect that only the first term in the summation on 
the right hand side of equation (IA6I) would need to be retained. If we define ipi as the angle between 
the gradient of the i-component (i.e., i = x,y,z) of Bo and e^, then when averaging over a large 
polarization map we have 



(B 0>i (x) [£ (e e ■ V) B 0jl (x)]> = £B 0>i (x) | VB ,i| (cos (</>*)) . 



(A7) 



But since is equally likely to be oriented in any direction over the whole map we have (cos {tpi)) = 
and the first order term of the Taylor expansion (i.e., equation |A7j ) cancels out. It therefore follows 
that the first non- vanishing term in the summation on the right hand side of equation (jA6|) is of 
second order with 



B (x) • 



oo 

£^-V) n B Q (x) 



n=l 



^(B (x)-(e £ -V) z B (x)^' 2 



(A8) 



when £ <C d. If we also assume stationarity for the turbulent magnetic field, then we define the 
autocorrelation of the turbulent field as 
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(BfB t (l)) = (Bt(x).B t (x + £)>, (A9) 
which, if we now define 5 as the correlation length for B t (x), has the following limits 

since the turbulent field is assumed uncorrelated over separations exceeding 5 and (B t ) = from 
the third of equations (IA2I) , Inserting equations (1A8|) and (|A9I) into equation (|A5I) we have 



(B (x) ■ B (x + 1)) ~ Bl (x) + ~ (b (x) • (e £ • V) 2 B (x)) t 2 + (B t • B t (£)} , (All) 
when f «i 

Using the assumed homogeneity in the fields' strength (i.e., equations [A3] ) the denominator 
of equation (|A4j) can be readily simplified to 



[<B 2 (x)><B 2 (x + £)>] 1/2 = <B 2 > 

= ( J B 2 + J B 2 + 2(B .B t )), 



which, with the fourth of equations (IA2|) . becomes 



[<5 2 (x)><B 2 (x + £)>] 1/2 = B 2 + <i? 2 >. (A12) 

If we further assume isotropy over space (i.e., A<1>3d (£) = A<3?3d {£)) and insert equations (|A11|) 
and (|A12p into equation (1A4|) we have 



[B?) - (Bt • B t (£)) - \ (Bp (x ) • (e, • V) 2 B (x) 



(cos [A$ 3 D W]) ^ 1 ~ 2 /K2V , (A13) 



when £ <^ d. For cases where A<3?3d is small equation (1A13|) simplifies to 



, , , NV 2 5 -(BfBt(Q) ( B (x) • (e, • V) B (x) 
A$1 D (*)) ~ LV t7 = - + » r^x , (A14) 
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still when t <C d. 

Examining equation (IA10I) we recover the behavior of the turbulent contribution to (A<J?| D (£)} 
(i.e., the first term on the right-hand side of equation |A14| ) described in Sj3] that goes from when 
t — > to a constant, which we now define as 6| D , when I > S. The data sets analyzed in this paper 
are such that t > 5 in all cases. We therefore find that the dispersion function is of the form 

<A*i D (^)>~6i D + m| D ^ 2 , (A15) 

with 

&3D " Bl + <5?> 

when 5 < £ <C d. Once again, we identify 63D with the constant contribution stemming from the 
turbulent field to the total angular dispersion, while the larger scale contribution due to variations 
in the large-scale field Bo is accounted for by the presence of a term proportional to £ 2 in equation 
([Al5jl . 

A. 2. Analysis in Two Dimensions 

The analysis presented above can still be used when we limit ourselves to two dimensions. 
This is needed in order to enable comparisons with polarization measurements, which only probe 
the plane-of-the-sky component, By, of the magnetic field. Defining e_|_ as the unit vector directed 
along the line-of-sight we have for the total magnetic field 

By =B-(B-e_L)e_L, (A16) 

and similar relations for Bo and B t . 

We need to evaluate, among others, the following autocorrelation 

(By (x) • By (x + £)) = (B (x) • B (x + £)) - ([B (x) • ej [B (x + £) ■ ej) , (A17) 

where the separation £ is now limited to the plane-of-the-sky. The last term on the right hand-side 
can be transformed to 



<[B (x) • e J [B (x + £) ■ ej) = <{[B (x) + B t (x)] • e ± } {[B (x + £) + B t (x + £)} • e ± }) 

= 5 2 ± + <[B t (x)-e_ L ][B t (x + £)-e ± ]>. (A18) 
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Using the same method that led to equation (|A12|) in the three-dimensional case we also have 



that 



BHx))(BHx + €)) =b1 1 + (b! 



1/2 



(A19) 



We now introduce the function 



(cos [A$ (£)}) 



(By (x) - By (x + l)) 
Bl (x)\ ( x + £) 



1/2 ' 



(A20) 



which upon inserting equations (lAlip . (IA17|) . (IA18j) . and (IA19I) with the condition of space isotropy 
becomes 



(cos [A$ {£)}) ~ 1 



B t ,|| • By (£)) - | (b (x) • (e, • V) 2 B (x)\ £ 



when I <C d and where 



(By • By {£)) = (B t • B t {£)) - ([B t (x) • ej [B t (x + £) • ej) 



If we further consider A<£> {£) to be small, then we find 



By - By {£)) B (x).(e r VrB (x) 



B l II + W| 



(A21) 



still when £ <C d and the displacement £ is limited to the plane-of-the-sky. 

For our data sets we have the further simplification that 5 < £ <C d and the dispersion function, 
equation (IA21I) . has then a form similar to its three-dimensional counterpart with 



(A$ 2 {£)) ~6 2 + m 



2,2 



(A22) 



where 



B l II + ( B l\ 



(A23) 
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is the quantity we evaluate through polarization measurements. Equation (|A23j) can be transformed 
to yield the ratio of the turbulent to large-scale magnetic field strength through 



1/2 



B, 



o,|| 



(A24) 



B. Data Analysis 

Data from the Hertz polari meter on the c l ouds studied here ha ve been previ ously published by 
Honrie et al.1 (hooi ) for OMC-1, iHonrie et all (l2002h for M17, and iKirhvl (£3) for DR21fMain ). 



Detail s on the instrum e nt as well as data acquisition and reduction can be found in iDowell et al 



t| 19981 ) and lKirby et al.l (|2005l ). respectively. The analy sis presented here i s per formed on a complet e 
re-reduction of the raw Hertz data using the method of lKirby et al.1 (|2005l ) and lDotson et al.1 ((2009). 
The data may differ slightly from that published in the references above. For our purposes we only 
include data which satisfy the P > 3ap criterion, where P is the polarization fraction and up its 
uncertainty. 

In each of the three objects the angle differences between each and every pair of data points 
are calculated as 



A$ 



Kl 



(Bl) 



and the corresponding distance between each point 



X; 



(B2) 



Note that lij = £ji so that a map with N data points contains only N(N — l)/2 distinct differences. 
Also note that |A<J>jj| is constrained to be in the range [0,90] degrees. 

These data are divided into separate distance bins with sizes corresponding to integer multiples 
of a single Hertz pixel-to-pixel separation (17'.'8); the first bin covers i\ < £ij < £2 (where 4 
corresponds to k pixels). Within each bin k we calculate the dispersion as the root-mean-square of 
the angle difference 



2 a/a 



for all lj~ < 



■v 



< 



(B3) 



The dispersion is corrected for measurement uncertainty within each bin according to equation 
(|3|). The uncertainty on each A^- follows from simply propagating the measurement uncertainties 
on both <3?i and $j through equation (|B1I) . The root-mean-square measurement uncertainties within 
each bin k are then given by 
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1 /2 

^M,fc = <(7 2 (A$i i )) fc / for all 4 < % < 4+1- 

The corrected dispersions are those plotted for the different bins in Figure [2j The error bars 
in Figure [2] are determined by propagating the measurement uncertainties for <5j and &j through 
equations (|B1I) and (lB3j) . most of these are too small to be seen in the figure, especially at the 
smallest displacements. 

For each object, the data are fitted to a linear model of the corrected square dispersion with 
respect to the square of the distance according to equation ([2|). In the discrete notation introduced 
in this section, the model is given by 

(A^) k -a 2 Mjk = b 2 + m 2 £l 

In order to ensure we are in the linear regime, the fits are limited to only the smallest three 
distance bins. Taking into account the uncertainties on the /a<E>?-\ , the least-squares solutions 
for the parameter b are given in Table [TJ 

REFERENCES 

Chandrasekhar, S., and Fermi, E. 1953, ApJ, 118, 113 

Crutcher, R. M., Troland, T. H., Lazareff, B., Paubert, G., and Kazes, I. 1999, ApJ, L121 
Crutcher, R. M., Nutter, D. J., and Ward-Thompson 2004, ApJ, 600, 279 
Curran, R. L., and Chrysostomou, A. 2007, MNRAS, 382, 699 

Dotson, J. L., Davidson, J. A., Dowell, C. D., Hildebrand, R. H., Kirby, L., and Vaillancourt, J. E. 
2009, submitted to ApJS 

Dowell, C. D., Hildebrand, R. H., Schleuning, D. A., Vaillancourt, J. E., Dotson, J. L., Novak, G., 
Renbarger, T., and Houde, M. 1998, ApJ, 504, 588 

Evans, N. J. II 1999, ARA&A, 37, 311 

Falceta-Gongalves, D., Lazarian, A., and Kowal, G. 2008, ApJ, 679, 537 

Frisch, U. 1995, Turbulence: The Legacy of A. N. Kolmogorov (Cambridge: Cambridge University 
Press) 

Girart, J. M., Rao, R., and Marrone, D. P. 2006, Science, 313, 812 



- 15 - 



Hoang, T., and Lazarian A. 2008, MNRAS, 388, 117 

Heitsch, F., Zweibel, E. G., Mac Low, M-M., Li, P., and Norman, M. L. 2001, ApJ, 561, 800 
Hiltner, W. A. 1949, ApJ, 114, 241 

Houde, M., Peng, R., Phillips, T. G., Bastien, P., and Yoshida, H. 2000, ApJ, 537, 245 

Houde, M., Bastien, P., Dotson, J. L., Dowell, C. D., Hildebrand, R. H., Peng, R., Phillips, T. G., 
Vaillancourt, J. E., and Yoshida, H. 2002, ApJ, 569, 803 

Houde, M. 2004, ApJ, 616, Llll 

Houde, M., Dowell, C. D., Hildebrand, R. H., Dotson, J. L., Vaillancourt, J. E., Phillips, T. G., 
Peng, R., and Bastien, P. 2004, ApJ, 604, 717 

Kleiner S. C, and Dickman, R. L. 1984, ApJ, 286, 255 

Kleiner S. C, and Dickman, R. L. 1985, ApJ, 295, 466 

Kleiner S. C, and Dickman, R. L. 1987, ApJ, 312, 837 

Kirby 2009, ApJ, in press 

Kirby, L., Davidson, J. A., Dotson, J. L., Dowell, C. D., and Hildebrand, R. H. 2005, PASP, 117, 
991 

Kudoh, T., and Basu, S. 2003, ApJ, 595, 842 

Lai, S.-P., Crutcher, R. M., Girart, J. M., and Rao, R. 2002, ApJ, 561, 864 

Lai, S.-P., Crutcher, R. M., Girart, J. M., and Rao, R. 2002, ApJ, 566, 925 

Lai, S.-P., Girart, J. M., and Crutcher, R. M. 2003, ApJ, 598, 392 

Lazarian, A., Vishniac, E. T., and Cho, J. 2004, ApJ, 603, 180 

Li, H., and Houde, M. 2008, ApJ, 677, 1151 

Myers, P. C, and Goodman, A. A. 1991, ApJ, 373, 509 

Ostriker, E. C, Stone, J. M., and Gammie, C. F. 2001, ApJ, 546, 980 

Padoan, P., Goodman, A. A., Draine, B. T., Juvela, M., Nordland, A., and Rognvaldsson, O. E. 
2001, ApJ, 559, 1005 

Schleuning, D. A. 1998, ApJ, 493, 811 

Vaillancourt, J. E. 2002, ApJS, 142, 53 



16 

Wiebe, D. S., and Watson, W. D. 2004, ApJ, 615, 314 



This preprint was prepared with the AAS IATgjX macros v5.2. 



-17- 




Fig. 1. — Dispersion: Idealized plots of the angular dispersion function, (A<I> 2 (^)) , between pairs 
of magnetic field vectors separated by displacements £, for values of £ <C d, with d the typical length 
scale for variations in the large-scale magnetic field (see Sj3j). Curve A: no measurement uncertainty; 
no turbulence. Curve B: with measurement uncertainty, <7m- Curve C: with turbulence. Curves D 
and E: accounting for correlation in polarization angles at displacements I smaller than the larger 
of the telescope beam (1.22X/D) (Curve D) or the turbulent correlation length 5 (Curve E). 
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Fig. 2.— Angular dispersion function, (A$ 2 (f)} 1/2 , for M17, DR21(Main), and OMC-1. The 
turbulent contribution to the total angular dispersion is determined by the zero intercept of the fit 
to the data at i = 0. The measurement uncertainties were removed prior to operating the fits to 
the corresponding data sets. The results are given in Table [TJ 



- 19 - 



Table 1. Results for the dispersion, the turbulent to mean magnetic field strength ratio, the line 

widths, and the mean field strength. 



Object 


6 a 


<^ 2 > 


1/2 /R b 
/-DO 


a(v) 


B C 




(deg) 






(km s _1 ) 


(mG) 


OMC-1 


8.3 ±0.3 


0.10 


±0.01 


1.85 


3.8 


M17 


10.4 ±0.6 


0.13 


±0.01 


1.66 


2.9 


DR21(Main) 


6.8 ± 1.3 


0.08 


±0.02 


4.09 


10.6 



a Turbulent dispersion (i.e., the dispersion limit as t — > 0). 
b Calculated with equation j7]). 

c Calculated with equation (jSJ), assumes a density of 10 5 
cm -3 and a mean molecular weight of 2.3. These estimates 
are not precise to better than a factor of a few. The process of 
signal integration through the thickness of the cloud and across 
the telescope beam inherent to the polarization measurements 
has also not been taken into account. 



